Instability of the Mott or Lieb-Wu insulator caused by an infinitesimal perturbation 
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The half-filled ground state of the Hubbard model in one dimension is studied by Kondo-lattice 
theory. Because of the Kondo effect, any insulating ground state with a complete gap open is 
unstable in the presence of an infinitesimal perturbation. This fact casts doubt on the claim by E. 
H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968) that the half-filled ground state is the 
Mott insulator. Though the claim is based on a rigorous result given by the Bethe-ansatz solution, 
the rigorous result is simply a necessary condition for the ground state being an insulator. 

PACS numbers: 71.10.Fd, 71.30.+h, 75.10.Kt 



I. INTRODUCTION 

Mott localization is still a crucial and present-day 
issue.— Whether magnetism is itinerant-electron or local- 
moment magnetism is directly related to physics of Mott 
localization. High-temperature (high- Tc) superconduc- 
tivity occurs in the vicinity of a phase of a local-moment 
type of antiferromagnetic (AF) insulator)^ Then it is an- 
ticipated that the mechanism of high-Tc superconductiv- 
ity is closely related to physics of Mott localization^"— 

The Hubbard model is one of the simplest effective 
Hamiltonians to study Mott localization. According to 
Hubbard's theory,^' - when the on-site repulsion U is large 
enough such that U > 2z|t|, where t is the transfer in- 
tegral between nearest neighbors and z is their coordi- 
nation number, the band splits into two subbands; the 
subbands are called the upper and lower Hubbard bands 
and a gap between them is called the Hubbard gap. Ac- 
cording to Gutzwiller's theory,—"— when the on-site U is 
large enough and electrons are almost half filled such that 
U > 2z\t\ and N c:i L, where N and L are the numbers of 
electrons and unit cells, respectively, electrons are renor- 
malized into heavy quasi-particles. A mid-band appears 
at the chemical potential between the upper and lower 
Hubbard bandsfli the mid-band is called the Gutzwiller 
band in this paper. According to Brinkman and Rice's 
theoryfiS when N — L a metal-insulator (MI) transition 
occurs at Uc — 2z\t\. In a metallic phase at U < Uc 
and U ~ Uc, the density of states (DOS) is of a three- 
peak structure, the Gutzwiller band between the upper 
and lower Hubbard bands. In an insulating phase at 
U > Uc, the Gutzwiller band disappears and the Hub- 
bard gap is a complete gap. According to these theories, 
the half-filled GS can be a paramagnetic type of Mott 
insulator, which is simply called a Mott insulator in this 
paper, only if U is large enough such that U > Uc- How- 
ever, any non-half-filled GS is a metal for any U. 

The GS can also be an insulator if an AF gap opens.i2 
In a high-T phase at T > Ep/kB, where is the ef- 
fective Fermi energy and fee is the Boltzmann constant, 
electrons behave like local moments even if no gap opens. 
If J7 > 2z\t\ and Ep < z\t\, the high-T phase is in a 
sense the Mott insulator. Then a local-moment type of 
AF insulator, in which the Neel temperature Tn is as 



high as Tn ^ E^/k^, is an AF type of Mott insula- 
tor. An itinerant-electron type of AF insulator, in which 
Tn ^ Ep/kQ, is not any type of Mott insulator. 

In spite of a lot of studies, the nature of the Mott in- 
sulator at T = K is not clear in the least except for the 
opening of the Hubbard gap. For example, it is contro- 
versial whether the MI transition at Uc is a first-order 
or second-order one, or a discontinuous or continuous 
one. According to Brinkman and Rice's theoryji^ either 
of the specific-heat coefficient and the spin susceptibility 
diverges as U Uc — 0- This suggests that the transition 
is continuous. However, no symmetry seems to be broken 
in the Mott insulator at U > Uc, residual entropy seems 
to remain and the third law of thermodynamics seems 
to be broken in it. This suggests that the transition is 
discontinuous. On the other hand, it is unlikely that the 
third law is broken when U is finite. Since the appear- 
ance of rigidity due to symmetry breaking can make the 
third law valid in an ordered state^l^ it is possible that 
a hidden order exists in the Mott insulator. It is con- 
troversial whether or not residual entropy remains in the 
Mott insulator. Under the present circumstances, it is 
desirable to critically examine whether the GS can be 
the Mott insulator, from the very beginning. 

When U is infinite and N = L in the canonical ensem- 
ble, each unit cell is occupied by a single electron with 
an arbitrary spin and there is no empty nor double occu- 
pancy. No electron is itinerant; no bosonic charge excita- 
tion is possible or the charge gap is infinitely large, while 
all the bosonic spin excitations are degenerate at the zero 
energy. This is a prototypic Mott insulator. Since its 
residual entropy is k^ In 2 per unit cell, even the GS of it 
is in a sense a high-T phase. When electrons are removed 
from it, empty occupancies or holes are itinerant Since 
electrons become more itinerant as U becomes smaller, 
any non-half filled GS must be a metal for any finite or 
infinite U, at least if no symmetry is broken. 

The Hubbard model in one dimension is particular: No 
symmetry can be broken in itf^ and the Bethe-ansatz 
solution for it was given by Lieb and Wui^i Denote the 
GS energy by Eq(N) as a function of TV. Two chemical 
potentials are defined for the GS of electrons: 



fi+iN)^EGiN+l)-EGiN), 



(1.1a) 
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for the addition of an electron and 

^I^{N) = EGiN)-EG{N^l), (1.1b) 

for the removal of an electron. Then a gap is defined by 

e,{N)=^i+{N)-ti^{N). (1.2) 

When N ^ L, eg{N) = 0. Then, any non-half-filled GS 
is a metalfi^ i.e., a Tomonaga-Luttinger (TL) electron 
liquidji^Ti^ When N = L, on the other hand, 

eg(L) > 0, (1.3) 

for U/\t\ > 0. Based on this rigorous result, Lieb and Wu 
claimed that an MI transition occurs at ?7c/|t| — and 
the half-filled GS is the Mott insulator for any nonzero 
U; this insulator is also called a Lieb-Wu insulator in this 
paper. 

When U is finite, the third law of thermodynamics is 
not broken in the Bethe-ansatz solution. It sounds cu- 
rious that the nature of the Mott or Lieb-Wu insulator 
is different from that of the prototypic Mott insulator, 
in which the third law is broken. The rigorous result 
(|1.3p is no sufficient condition for the GS being an insu- 
lator but simply a necessary condition for it; e.g., eg{L) 
is nonzero in a metallic fine particle because of the long- 
range Coulomb interaction. It is desirable to reexamine 
whether the half-filled GS in one dimension is really an 
insulator, in particular, in the grand canonical ensemble; 
electrodes are necessary to measure the static conduc- 
tivity of a system and the electrodes play a role of an 
electron reservoir for the system. 

Not only Eq{L) but also eg(i) as a function of U arc 
singular aiU = ;^^i^^ they cannot be expanded in terms 
of U. Thus, eg(i) is extremely or singularly small for 
U/\t\ < 1: 

^lim^eg(i)/(C//K|)'' = 0, (1.4) 

for any real 77 > 0, though it is as large as the gap given 
by Hubbard's theory, or the Hubbard gap, for U > 2\t\. 
Since no symmetry is broken, no rigidity can appear in 
the gapped phase characterized by eg(i) > 0. Since it 
is anticipated that such a non-rigid gapped phase is un- 
stable at least in the limit of U /\t\ — >■ in the presence 
of a perturbation such as one due to an electron reser- 
voir, it is interesting to examine how large Uc can be as 
a function of the strength of the perturbation. 

In the grand canonical ensemble, the averaged electron 
number {N) can be an irrational one even in the limit of 
T — > OK; thus, it is a continuous function of the chemi- 
cal potential. This means that quantum fluctuations of N 
are effectively considered in the conventional formulation, 
though an electron reservoir is only implicitly considered 
in it. Quantum fluctuations of N are effective at least 
when (iV) is a non-integer, e.g., even when {N) = L± 0. 
However, they may or may not be effective when exactly 
{N) = L: If they are not effective, the half-fllled GS is 



composed of only components oi N = L electrons, so 
that it is anticipated that the gapped phase is rigid and 
Uc/\t\ — 0, i.e., the half-fllled GS is an insulator for any 
nonzero U in the grand canonical ensemble; though the 
GS is a metal even for {N) = L ± 0, it is an insulator for 
exactly {N) = L?^ On the other hand, if quantum fluctu- 
ations of N are still effective even when exactly {N) = L, 
the half-fllled GS is composed of not only components of 
N = L electrons but also those oi N ^ L electrons. Since 
any non-half-fllled GS oi N ^ L electrons is a metal in 
the canonical ensemble, it is anticipated that the gapped 
phase is not rigid and f/c/|^| — +00, i.e., the half-fllled 
GS is a metal for any flnite U in the presence of an elec- 
tron reservoir or in the grand canonical ensemble. It is a 
crucial issue which is the truth. 

Based on Kondo-lattice theory (KLT);^"— it was 
shown in a previous paper^^ that the GS cannot be the 
Mott insulator for any fllling {N) in the presence of a 
simplifled electron reservoir, though it is never any ac- 
tual type of reservoir or electrodes. In this theory^^l the 
Kondo effect plays a crucial role in causing any insulat- 
ing GS that is possible within the constrained Hilbcrt 
subspace where no symmetry is allowed to be broken to 
be unstable. The main purpose of this paper is to show 
that because of the Kondo effect the half-fllled GS in one 
dimension cannot be the Mott or Lieb-Wu insulator in 
the presence of an inflnitesimal perturbation, though it 
is a phenomenological one that guarantees the average 
{N) to be able to be an irrational number. Another pur- 
pose is to examine or discuss whether the half-fllled GS 
is the Mott or Lieb-Wu insulator in the absence of the 
inflnitesimal perturbation. 

This paper is organized as follows: Preliminary is given 
in Secim Result is given in Sec. lIIII Discussion is given 
in Sec. lIVI Conclusion is given in Sec.|Vl In Appendix 
El it is shown that the resonating-valence-bond (RVB) 
mechanism^* stabilizes a metallic GS rather than the 
Mott insulator. In Appendix|Bl it is shown that the 
Kondo effect is irrelevant in a non-lattice model or a con- 
tinuous model. In AppendixICl it is shown that the Mott 
insulator in inflnite dimensions is also unstable in the 
presence of the inflnitesimal perturbation. 

II. PRELIMINARY 
A. Kondo-lattice theory 

We consider the Hubbard model with a flnite onsite 
repulsion U in one dimension: 

n^-tJ2 dldja + UJ2 '^^t"^. (2-1) 

{ij)a i 

where riia = d\^dic and (ij) stands for nearest neighbors; 
the other notations are conventional. An inflnitesimal 
perturbation is phenomenologically considered. 

Since no symmetry can be broken in one dimension)^ 
the Hilbert space is constrained within the subspace 
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where no symmetry is allowed to be broken. Then the 
thermal Green function of electrons is given by 

G.(z£/,fc) = i (2.2) 

tSi+n-E{k)~T,^ {lei ,k) + i5- — r 

where ei = (21 + l)7rfcBT, with I being an integer, is a 
fcrmionic energy, /i is the chemical potential. 



E{k) = -2icos(fca), 



(2.3) 



with a being the lattice constant, and S^(iei,fc) is the 
self-energy. The term of i5ei/\ei\ is due to the phe- 
nomenological perturbation; the limit of (5 ^ or (5 = 0+ 
is assumed in this paper. The main purpose of this paper 
is to show that the half-filled GS cannot be an insulator 
in the presence of 5 = 0+ . 

In general, the self-energy is decomposed into the 
single-site Scr(«ei) and the multisite AEo-(«ez,fc): 

^a{tei,k) = t^iiei) + AJ:„{iei,k). (2.4) 

According to the Kondo-lattice theory (KLT)^^"^^ or dy- 
namical mean- field theory (DMFT)^ calculating the 
single-site T,^{iei) is mapped to a problem of self- 
consistently determining the Anderson model (AM) and 
solving it, i.e., a self-consistent Kondo problem. 
We consider the AM defined by 



(2.5) 



where fida = d]jda] the notations here are also conven- 
tional except for the use of tildes. The parameters of the 
AM are determined in a way such that the self-energy 
of it is equal to the single-site Yia{i£i) of the Hubbard 
model. Then the Green function of the AM is given by 
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1 r+°° A(e) 
^£^+/2-ed-Scr(^e^) / de -. 



-, (2.6) 



where jl is the chemical potential and 

K{e) = lY,\V^\H[e + fi-E,{\,)\. (2.7) 
^ k 

The Fermi surface (FS) is defined by // = i?c(k). If the 
FS exists, then A(0) > unless |Vkp vanishes on the 
whole of the FS. 

For the Hubbard model, a single-site Feynman diagram 
is defined as one such that it is only composed of on- 
site interaction and site-diagonal electron lines in the site 
representation: U and 



Ra{i£i) = ;^ X! G^{iei,k). 

k 



(2., 



For the AM, there is a Feynman diagram corresponds 
to the single-site one. It is composed of on-site interac- 
tion and electron lines: U and Ga{i£i)- The condition to 
determine the AM is so simple that U — U and 



Ra{iei) = G„{iei). 



(2.9) 



The condition (|2.9|) is equivalent to a set of /i = /i — 
and 



A(£) = Im t„{e + iO) + l/R„[e + iQ) 



(2.10) 



where ^^(e + iO) and Ra{e + iO) are the analytical con- 
tinuations of Y,„{iei) and Ra{iei), respectively, from the 
upper half plane onto the real axis. Since U — U and 
fi — fl — id are trivial, Eq. (|2.10|) is a practical mapping 
condition to determine the AM. 

If the multisite AY,cr{i£i, k) is ignored in this theoreti- 
cal framework, the theory is simply within the supreme 
single-site approximation (S'^A),'^^"'^'^ which is rigorous 
in infinite dimensions within the constrained Hilbert sub- 
space where no symmetry is allowed to be broken. Either 
the DMFT or dynamical coherent potential approxima- 
tion (DCPA)M is also within the S^A. In the KLT of 
this paper, a set of A(e), ^^(e -f iO), and AI]cr(£ + iO, k) 
should be self-consistently considered and determined to 
satisfy Eq. ((2im . 

Though the contribution of a Feynman diagram for 
the self-energy may diverge as T OK and |e| or 
|£;| — >■ for finite I, it is finite for T > OK because no 
symmetry is broken in one dimension. Therefore, it is 
certain that the mapping is possible at least for T > K. 
When the summation of all the Feynman diagrams is 
followed by taking the limit of T K, the mapping is 
possible even in the limit of T ^ OK. 

Because of Eq. ([2^ . the DOS of the Hubbard model 
is equal to that of the mapped AM: 



p(£) = — -y^lmG„(e + iO,k) (2.11a) 

k 

^ ~ (2.11b) 



= IniG^(£ + iO). 

TT 



If /9(0) > 0, the GS is a metal. If p(0) = 0, the GS is a 
non-metal: an insulator or a zero-gap semiconductor. 
The averaged electron number is given by 



where 



{N) = I def{e)p{e), 



stands for the statistical average, and 

1 



/(e) = 



(2.12) 



(2.13) 



(2.14) 
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In the thermodynamic hmit of L — ^ +00 followed by the 
limit of ^ — ?> 0, {Af) can be an irrational number and is a 
continuous function of fj, even in the limit of T ^ K. 

When fj, = U/2, both the Hubbard model and the 
mapped AM are half-filled and symmetrical: 

^1 = jl-id = U/2 = U/2, (2.15) 

{M)/L = (nt + = 1, p{s) = p{-e), 

A(e) = A(-e), (2.16) 

and so on. It follows from Ea. (|2.16p that 

/ TT^ =«A(0). (2.17) 

In a self-consistent solution, I](j(e -I- iO) or I](j(e -I- iO, k) 
may be continuous or discontinuous at e = 0. In either 
case, 

jl- id-Ref^ai+iO) ^0, (2.18) 

and 

fi- E{±kF) -ReT.„{+iO,±kF) =0, (2.19) 

where fcp = 7r/2a, so that £'(±fcF) = 0. 

In this paper, the half-filled GS in the symmetrical 
Hubbard model is mainly studied; < U/\t\ < +00, 
L +00, and T OK are assumed. 

B. Kondo effect 

The Kondo problem for the s-d model and the AM has 
already been solved)^"— The single energy scale of fceTk 
only appears in low-energy phenomena; Tk is called the 
Kondo temperature. In this paper, it is defined by 

fcBrK = l/[Xs(0)]y^OK> (2-20) 

where Xs (0) is the static susceptibility of localized elec- 
trons in the AM; in this paper, the susceptibility is de- 
fined in a way such that it does not include the conven- 
tional factor (7^^1/4, where g and /ie are the g factor 
and the Bohr magneton, respectively. The Kondo tem- 
perature depends on A(e) or A(0), which is defined by 
Eq. ([271): If the FS exists and 

A(0) > 0, (2.21) 

then Tk > OK. If A(0) = 0, 7k > OK or Tk = OK. The 
nature of the GS depends on Tk: If 7k > OK, the GS is 
a singlet and a normal Fermi liquid (FL). If Tk = OK, 
the GS is neither a singlet nor a normal FL; residual en- 
tropy remains and localized electrons behave like a local 
moment at any temperature T. 

According to the FL theory by Yamada and Yosida^i^ 
and the Bethe-ansatz solutior^ir— for the AM, the adia- 
batic continuation as a function of U holds in the region 



of < C7/[7rA(0)] < +00; i.e., the transition between the 
weak-coupling region of {//[7rA(0)] ^ 1 and the strong- 
coupling region of ?7/[7rA(0)] ^ 1 is only a crossover 
and there is no finite critical U. In this paper, therefore, 
physics in the AM for any U is simply called the Kondo 
effect. 

In the self-consistent Kondo problem of this paper, 
eventual A(e) and Tk are renormalized by intersite ef- 
fects, so that they should be self-consistently determined 
with the multisite AT^^{e + iO, fc) to satisfy Eq. (|2.10l) . 
The nature of the GS of the mapped AM depends on the 
eventual and self-consistent A(0) and Tk- 



III. RESULT 

A. Case of 5 = 0+ 

According to Eq. (|2.10p . it follows that 

A(e) > S, (3.1) 

for any e; the proof of Eq. p.ip is in parallel with that 
in Appendix of the previous paper. It should be noted 
that this inequality holds even if the self-energy S(j(e-|-iO) 
or AY,cr{s + iO, k) is anomalous. If (5 > 0, then A(e) > 0; 
thus, Tk > K for any self-consistent solution, though 
Tk may be infinitesimally low for S — O'^ . 

When i5 = 0+, the GS of the mapped AM is a normal 
FL because of the Kondo effect. In the presence of an 
infinitesimal magnetic field H in the AM, 'E^{e + iO) of 
the normal FL in the AM can be expanded in terms of 
e/(fcBTK) in a normal way such thali^Si^ 

i:^(e + iO) = ^U+{l- 4)^)e + (l - 0,) ^ag^iBH 

- i4>r£VikBTK) + O (e^) , (3.2) 

where <^^, 4>s, and 4>t are positive constants. The Wilson 
ratio is defined by 

= 4>s/4>j, (3.3) 

which is used in AppendixfAl In this section, iJ = is 
assumed. It follows that 

Gaie + iO, k) = l/[0-,e - E{k) + i^re^/{kBTK) 

- AS^(£ + i0,A:) + O(e3)]. (3.4) 

In one dimension, the charge-spin separation occurs 
in pair or bosonic excitationsi^Ti^i and a single-particle 
fermionic excitation is in a sense a composite of charge 
and spin excitations j^^'^^ however, it is doubtful or at 
least unclear whether charge-spin separation occurs in 
the fermionic single-particle excitation. In any case, it is 
inevitable that the multi-site AScr(e + iO, k) is anomalous 
at e = in one dimension. 
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If AEo.(e + iO, k) is continuous and finite at e = 0, it 
follows from Eqs. (l2.11ap . ([^1^ . and ([331) that 



and 



lim p{e) > 0, 

e— >-0 



p(0) > 0, 



(3.5a) 



(3.5b) 



for the Hubbard model. It also follows from Egs. (|2.6p . 
([2llbl) . (l2Tfl) . ([2J81) . and (l33|) that 



p(0) - l/[7rA(0)] > 0, 



(3.6) 



for the AM. Since /cbTk is the only low-energy scale in 
the AM, Eq.dnH]) implies that p{e) > for at least \e\ < 
k^TK- Thus, the GS is a metal in this case. 

Because of the Kramers-Kronig relation, if one of the 
real and imaginary parts of AI]o-(e + iO, k) is discontin- 
uous at e = the other is divergent at e = 0. Even if 
AI]cr(£ -I- zO,fc) is divergent at e = 0, no complete gap 
opens because of the imaginary term of i^^-e^/ (fceTx) in 
the single-site self-energy but only a zero-gap can open. 
For example, if AI]cr(£ + *0, k) has a pole at £ = — iO in 
a way such that 



AS^(£ + iO,fc) = ci/(£ + iO) 



then 



p(£)- [4>r/{c\k^T^)\t 



O 



(3.7) 



(3.8) 



i.e., only a zero-gap opens. Even if AJ^a-{£ + *0, fc) is di- 
vergent at £ = 0, the GS cannot be an insulator with a 
complete gap open but it can be only a zero-gap semi- 
conductor. 

When U is finite and 5^0+, the half-filled GS is 
a metal or a zero-gap semiconductor and it can never 
be any insulator with a complete gap open, so that 
[/c/l^l = +00. The only approximation used to obtain 
this conclusion is that the expansion (|3.2I) is made use 
of. If the expansion coefficients of (f>s, and (jj-r are 
rigorous, the expansion becomes rigorous in the limit of 
£/(fcB'7k) — i> 0. Essentially, no approximation is used 
to obtain the conclusion; thus, it is unquestionable and 
absolute. 



B. Case of absolutely 5 — 

When U is finite and 5 — 0, there are two possibil- 
ities for Tk: Tk > ok and Tk = OK. Provided that 
Tk > OK, the expansion p.2p is relevant, so that the 
single-site self-energy is normal. As proved in Sec. lIII Al 
no complete gap opens because of the Kondo effect; the 
GS is a metal if the multisite self-energy ASo-(£ -I- iO, k) 
is continuous and finite at £ = or a zero-gap semicon- 
ductor if AEo-(£ + *0, k) is discontinuous or divergent at 
£ = 0. 



Provided that Tk = K, the expansion (13.21) is irrele- 
vant, so that not only the multisite self-energy but also 
the single-site self-energy may be anomalous. Then, the 
total self-energy Yia{s+iQ, k) may be so anomalous that it 
has a pole just below the real axis in the vicinity of £ = 
and its imaginary part ln\Yi„{e + iO, k) is absolutely zero 
for at least —eQ/2 < e < €q/2 except for the existence 
of the delta function at the singular point where the pole 
exists; e.g.. 



I]^(£ + iO,fc) = 



ci(fc) 



£ + iO TT 



1 A(e;.) 



(3.9) 



where A{e]k) = for at least -cq/Z < e < cg/S. A 
complete gap as large as eq can only open in such an 
anomalous case. Thus, the GS can be an insulator only 
if Tk = K; Tk = K is a necessary condition for the 
GS being an insulator with a complete gap open.^^ Essen- 
tially, no approximation is used to obtain this conclusion 
either; thus, it is also unquestionable and absolute. 

The only possible insulating GS is one such that it 
is characterized by Tk = K. What is proved in Sec. 
nil Al is simply that the insulating GS is unstable if the 
infinitesimal 5 = 0+ is once introduced. 



IV. DISCUSSION 

It is proved in Sec. lIIII that the half-filled GS is no 
insulator in the presence of the phenomenological 5 = . 
In this section, it is examined or discussed from various 
points of view whether or not the half-filled GS is an 
insulator in the absence of <5 = 0+ . 

In field theory based on the grand canonical ensem- 
ble, the analytical continuation from the upper half-plane 
onto the real axis is made in a way such that 



(4.1) 



with (5ac — 0^. This S^c = 0^ plays a similar role to 
(5 = O'^ in this paper. If a simplified electron reservoir is 
considered as in the previous paper, the thermal Green 
function is given instead of Eq. (|2.2p by 



Gaiiei, k) 



1 



iei +ii — E{k) - Y.cr{iei, k)+r{iei 



(4.2) 



where T(iei) is due to the hybridization with the reser- 
voir. If the FS exists in the reservoir, it follows that 



Imr(+iO) > 0. 



(4.3) 



This Imr(-fiO) > also plays a similar role to 5 = 0+ in 
this paper. It is plausible that Tk > K and the third 
law of thermodynamics is not broken even in the absence 
of (5 = 0+, unless U is infinite. Thus, it is desirable to 
critically examine whether the half-filled GS is really an 
insulator with a complete gap open in the absence of the 
infinitesimal phenomenological (5 = 0"*" 
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Not only the result of Sec. lIIII but also the argument 
above in this section cast doubt on whether Eq. (|1.3p or 
eg(L) > is really a proof that the half- filled GS for 
(5 = is an insulator. When Tk — OK, the residual 
entropy of localized electrons in the AM is ©(fee), so that 
the single-site residual entropy of the Hubbard model is 
also 0{k^) per unit cell. If an insulating GS is possible 
when 5 = 0, it should be one such that Tk = K and 
the third law of thermodynamics is broken in it. On 
the other hand, the third law is not broken in the Bethe 
ansatz solution, which implies that Tk > K even for 
(5 = 0. Thus, it is doubtful whether the half-filled GS of 
the Bethe-ansatz solution is the only possible insulating 
GS, which is characterized by Tk = OK. The rigorous 
result of eg(T) > is no proof that the half- filled GS of 
the Bethe-ansatz solution is the Mott insulator. 

When an electron or hole is added to the Hubbard 
model in the canonical ensemble, the whole of it remains 
within the Hubbard model. When it is added to the 
Hubbard model in the grand canonical ensemble and if 
quantum fluctuations of M are at least effectively consid- 
ered in a way such that (A/") can be an irrational number, 
not only a process where the whole or almost whole of it 
remains within the Hubbard model but also a process is 
possible such that only an infinitesimally small fraction 
of it remains within the Hubbard model and the almost 
whole of it escapes to a reservoir. The single-particle ex- 
citation spectrum p{£) in the grand canonical ensemble 
can be different from the spectrum of adding or removing 
an electron in the canonical ensemble. Thus, £g{L) > 
is no sufficient condition for the opening of a complete 
gap in p(e) but simply a necessary condition for it. 

The total lattice of the Hubbard model can be divided 
into sub-lattices in a way such that a central one is sur- 
rounded by the others. The others can play a role of 
an electron reservoir for the central sub-lattice. Denote 
the numbers of electrons and unit cells within the central 
sub-lattice by N' and L' , respectively, and the averaged 
number of N' by {N'): {N') = L. When U is finite, 
quantum fluctuations of N' develop to a greater or lesser 
degree and they never vanish.— The half-filled GS on the 
central sub-lattice is therefore composed of components 
of various N' such as N' = L' , L' zL 1, L' ± 2, and so on. 
Any non-half-filled GS, where N' ^ L\ is a metal. Since 
the half-filled GS on the central sub-lattice includes such 
non-half-filled metallic components, it is likely to be a 
metal. Since any sub-lattice can be a central one, the 
half-filled GS on it is also likely to be a metal. The half- 
filled GS on the total lattice is composed of such half- 
filled metals on the sub-lattices. Thus, the half-filled GS 
is likely to be a metal even in the canonical ensemble. 

According to a perturbative treatment,~ when 1 < 
C//|t| < -|-oo, the superexchange interaction arises from 
quantum fluctuations allowing empty and double occu- 
pancies. The same one can also be derived by a treat- 
ment of field theory] "^^1^° Since it arises from the virtual 
exchange of a pair excitation of an electron in the upper 
Hubbard band and a hole in the lower Hubbard band. 



it works in either a metallic or insulating phase. In the 
limit of U -l-oo, its constant is given by 

J = -U^jU, (4.4) 

between nearest neighbors; when U l\t\ > 1, Eq. (|4.4|) can 
be approximately used. The superexchange interaction 
causes not only the development of AF spin fluctuations 
but also the stabilization of a singlet GS by the forma- 
tion of an itinerant singlet or a resonating valence bond 
(RVB) on each pair of nearest neighbors, as studied in 
Appendix [X] This stabilization mechanism is simply the 
RVB mechanism proposed by Fazekas and Anderson for 
the Heisenberg model on the triangular lattice^^ The 
stabilization energy is 0(| J|) per unit cell; it is also the 
quenching energy of a spin at a unit cell. Then, 

fc(0)]T^0K=O(l/|J|). (4.5) 

From Eqs. (E^D]) and (03]), fceTK = 0{\J\). The fact that 
Tk > K for any finite U also casts doubt on whether 
the half-filled GS is an insulator when U is finite. 

In the Heisenberg limit oi U/\t\ — +oo with J = 
—At^/U kept constant, the half-filled Hubbard model is 
reduced to the Heisenberg model. The GS of the Heisen- 
berg model is a singlet or a doublet according to the 
number of spins, even or odd, unless J — Q. This is a 
proof or at least evidence that the GS degeneracy of the 
Hubbard model is never infinite even in Heisenberg limit. 
This also casts doubt on whether the half-filled GS is an 
insulator in the Heisenberg limit. 

It is obvious that Tk > 0, p(0) > 0, and A(0) > for 
the metallic GS while Tk = OK, p(0) = 0, and A(0) = 
for the possible insulating GS. If an MI transition oc- 
curs at Uc, the MI transition is inevitably a discontinu- 
ous one: Residual entropy jumps at Uc- Since Eq. (|3.6p 
is satisfied in the metallic GS, one or both of p(0) and 
A(0) also jump at Uc- It is shown in Appendix El that 
Uc/\t\ = +00 if the RVB term is only considered as the 
multisite self-energy beyond the S^A within the KLT; 
p(0) continuously vanishes as U/\t\ — > -l-oo while A(0) 
diverges as U/\t\ +oo and A(0) = for U/\t\ = +oo. 

Physical properties of a spin liquid in the Heisenberg 
model in one dimension can be mapped to those of a TL 
liquid, i.e., the GS of the Heisenberg model is a TL spin 
liquid^-^ This fact implies that the half-filled GS of 
the Hubbard model is not the Mott insulator but a TL 
electron liquid for any finite U even in the Heisenberg 
limit. When U/\t\ > 1 or /cbTk = 0(| J|), the TL electron 
liquid is a type of RVB electron liquid because it is mainly 
stabilized by the RVB mechanism. 

If no symmetry is broken in two dimensions and higher, 
e.g., because of frustration and low dimensionality or 
because of forced constraint of the Hilbert space, it is 
straightforward to extend the analyses in this paper to 
two dimensions and higher. When 5 = 0+, the GS can 
never be an insulator with a complete gap open even in 
the Heisenberg limit. If no symmetry is broken in a spin 
liquid in the Heisenberg model, it is probable that physi- 
cal properties of the spin liquid in the Heisenberg model 
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can also be mapped to those of an electron liquid in the 
Hubbard model, as in one dimension, except for conduc- 
tivity; e.g., those of the RVB spin liquid in the Heisen- 
berg model on the triangular lattice^^ must be mapped 
to those of an RVB electron liquid in the Hubbard model 
on the triangular lattice. Any type of spin liquid, in 
which the third law is not broken and the bosonic spin- 
excitation spectrum is normal, should be distinguished 
from the Mott insulator discussed in this paper, in which 
Tk = K and the third law is broken; Tk = K means 
or at least implies that there is degeneracy at the zero 
energy in the bosonic spin excitation. 

The Hubbard model in one dimension becomes the 
Tomonaga-Luttinger (TL) model, when E(k) given by 
Eq. (|2.3|) is replaced by the dispersion relation composed 
of two linearized branches defined by 

EThik) = Eikp) ± 2|i|asin(fcFa)(A: =F kp)- (4.6) 

The TL model can be mapped to a boson model; the 
mapping is rigorous only if the cutoff for Eq. (|4.6p is in- 
finitely large. In the boson model, the charge and spin 
parts are separable from each other, i.e., the charge-spin 
separation occurs J^^—i^ When kp = ■jT/2a, a gap opens 
in the charge excitation because of umklapp processes. 
The charge gap is as large as 

Ap = Ua/{2T:a), (4.7) 

for a particular U, where l/a is the cutoff discussed 
above; however, it is not clear whether Ap as a function 
of U is analytic or singular a.t U = and whether Ap 
corresponds to eg{L). Even in the gapped phase where 
the charge gap Ap opens, the spin-excitation spectrum is 
normal: No gap opens in it and no degeneracy exists at 
the zero energy. If Tk is a well-defined property in the 
boson model, this normal spectrum means that Tk > K 
for the gapped phase. In the possible insulating GS of 
the Hubbard model, on the other hand, Tk = K and, 
presumably, there exists degeneracy at the zero energy 
in the bosonic spin-excitation spectrum. The third law 
of thermodynamic is not broken in the gapped phase of 
the boson model, while it is broken in the possible insu- 
lating GS of the Hubbard model. Thus, the nature of the 
gapped phase of the boson model is different from that of 
the possible insulating GS of the Hubbard model, which 
casts doubt on whether the opening of the charge gap Ap 
in the boson model corresponds to the opening of a gap 
in the fermionic p{e) of the Hubbard model. 

The argument in the previous paragraph also casts 
doubt on whether Tk is a well-defined property in the 
boson model and the TL model. In the Hubbard model, 
the Kondo effect plays a crucial role in causing any insu- 
lating GS to be unstable, at least, when 6 — 0+. There 
are two possible scenarios for the relevance of the Kondo 
effect in the TL model: In a scenario, it is a relevant ef- 
fect, so that Tk and Scr(£-l-«0) are well-defined properties 
and no gap opens if ^ = 0^ , as in the Hubbard model. In 
the other scenario, it is not a relevant effect, so that nei- 
ther Tk nor 1^^(6 + 10) is any well-defined property and a 



gap opens even if 5 = 0+; the opening of the charge gap 
in the TL model does not correspond to the opening of 
a gap in the fermionic p{e) of the Hubbard model. The 
Hubbard model is a lattice model, while the TL model is 
essentially a continuous model. Since the Kondo effect, 
which is a single-site effect, is peculiar to a lattice model, 
the latter scenario is likelier than the former one; e.g., it 
is shown in Appendix|B] that the Kondo effect can play 
no role in a continuous model. The equivalence between 
the Hubbard model and the TL model, in particular, as 
regards the opening of the charge gap should be criti- 
cally reexamined from a point of view whether or not the 
Kondo effect is a relevant effect in the TL model. 



V. CONCLUSION 

The half-filled ground state of the Hubbard model 
is studied by the Kondo-lattice theory within the con- 
strained Hilbert subspace where no symmetry is allowed 
to be broken. The Kondo temperature Tk or fceTk is 
the energy scale of single-site quantum spin fiuctuations. 
Provided that Tk > OK, the single-site self-energy for 
electrons is normal, so that the ground state is never any 
insulator with a complete gap open; it is a metal if the 
multisite self-energy for electrons is continuous and finite 
at the chemical potential or a zero-gap semiconductor if 
it is discontinuous or divergent at the chemical poten- 
tial. Thus, Tk = K is a necessary condition for that 
the half-filled ground state is an insulator with a com- 
plete gap open. On the other hand, any ground state 
characterized by Tk = K is infinitely degenerate in the 
thermodynamic limit, unless the electron filling is empty 
or completely filled. Thus, the only possible insulating 
half-filled ground state is one such that Tk = K and 
the third law of thermodynamics is broken in it. 

In the presence of an infinitesimal phenomenological 
perturbation, definitely Tk > K, unless the onsite U is 
infinite. The only possible insulating half-filled ground 
state, which is characterized by Tk — OK, is unstable 
in the presence of the infinitesimal perturbation, though 
the prototypic Mott insulator for infinite U is stable. 

Though Lieb and Wu's claim that the half-filled ground 
state in one dimension is the Mott insulator for any 
nonzero U is based on the rigorous result given by the 
Bcthe ansatz solution, the result is not a proof that the 
half-filled ground state is an insulator with a complete 
gap open because the result is not a sufficient condi- 
tion for it but simply a necessary condition for it. Since 
the third law is not broken in the Bethe ansatz solu- 
tion, which implies that Tk > K even in the absence 
of the infinitesimal perturbation, it is doubtful whether 
the half-filled ground state is the only possible insulating 
ground state, which is characterized by Tk = K. It is 
doubtful whether it is the Mott insulator. 

The Hubbard model is a lattice model. Any single- 
site property of the Hubbard model can be mapped to 
its corresponding one in the Anderson model that is 
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self-consistently determined by the Kondo-lattice theory. 
When Tk > OK, the Kondo effect plays a crucial roles 
in causing the possible insulating ground state to be un- 
stable. Since the Kondo effect, which is a single-site ef- 
fect, is peculiar to a lattice model, it is unlikely that 
the Kondo effect is a relevant effect in the Tomonaga- 
Luttinger model, which is essentially a continuous model. 
A gap can only open in the fermionic single-particle ex- 
citation spectrum of the Hubbard model if Tk = K and 
the third law is broken, while a gap opens in the charge 
excitation spectrum of the Tomonaga-Luttinger model, 
or the boson model to which the Tomonaga-Luttinger 
model is mapped, even if the third law is not broken. 
Thus, it is doubtful whether the opening of a gap in the 
charge excitation spectrum of the Tomonaga-Luttinger 
model or the boson model is a proof that a complete gap 
opens in the fermionic single-particle excitation spectrum 
of the Hubbard model. It should be reexamined whether 
the mapping of the Hubbard model in one dimension to 
the Tomonaga-Luttinger model and the boson model is 
really relevant as regards the opening of the gap in the 
bosonic charge-excitation spectrum. 
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Appendix A: RVB stabilization effect 

In this Appendix, one dimension, {J\f) — L, U/\t\ ^ 
1, and (5 = 0+ are assumed and the RVB term is only 
considered as the multisite self-energy beyond the S"^A 
within the KLT. The RVB term is the Fock-type term of 
the superexchange interaction: 



AJ:r^H^e,,k) = ^J2 cr'^''^) 
X G|,';)(^e„,p)e-"°^ (Al) 



where 



Js{q) = 2Jcos{qa), 



(A2) 



with J given by Eq. (|4.4p , cr'^'^ is the aa' component of 
the Pauli matrix cr — (cr^, ay, cr^), and Glj\i£n,p) is the 
low-energy coherent part of the total Ga-{iSn-,p), which 
is given (|2.2|) . Since the energy dependence of the su- 
perexchange interaction is ignored, the coherent part is 
only considered in Eg. (|A1[) : thus, Eq. (jAll) is only valid 
for \ei\ < fceTk if and only if Tk > OK. According to 
the Ward identity,— 0^ defined by Eq. p.2p is the static 
part of the reducible single-site three-point vertex func- 
tion in spin channels; 0s is approximately used as its low- 
energy dynamical part. Since multisite or intersite effects 



are perturbatively considered based on the S^A;^— the 
two single-site vertex functions should be included in Eq. 
there is no double counting in this treatment. 
When ^ = 0+, the expansion p.2p is relevant. Then, 
the coherent part is given by 



1 



(A3) 



and the RVB term is given by 



AE';^^^Hiei,k) = -J4>^W^Ecos{ka). (A4) 
Here, Wg is the Wilson ratio defined by Eq. p.3p and 

S= i^cos(fca)/[e(fc)], (A5) 



where /(e) is defined by Eq. (|2.14p and 

^(fc) = -2{t*/4>^) cos(fca), (A6) 

where 



In the limit of T K, S = I/tt = 0.31831 ••• . 
The coherent part is simply given by 



(A7) 



Gi-Htsi,k) 



1 1 
iei - C(fc) ■ 



(A8) 



The bare dispersion relation E{k) is renormalized into 
^(fc) by the RVB mechanism. The DOS of the coherent 
part is given by 



P 



(A9) 



which is simply the DOS of the Gutzwiller band; 

p(^)(e) = (ap/|r|)[l + 0(£2)], (AlO) 

with Up = 0(1), for |e| < \t*\/^~f. The bandwidth of ^(A:) 
or p^'^^(e) is 

W* ^2\t*\/~^^, (All) 
which corresponds to /cbTk. It should be noted that 

W* ■i\J\W'^E/A, (A12) 

if <f)~f +00, i.e., W* is nonzero even if 0^ — > -|-(X), 
while p('^'(e)|i| — ^ if — -|-oo. It is interesting that, 
unless J = 0, the bandwidth W* of the Gutzwiller band 
is nonzero even though the band itself is vanishing. 

The expansion coefficients cf)^ and 0s , which are single- 
site properties, should be self-consistently calculated with 
the RVB term. In principle, this problem should be 
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solved by self-consistently determining and solving the 
AM, as discussed in Sec. Ill Al However, we take a qualita- 
tive approach to this problem. If the mapping condition 
is satisfied, the probabilities of empty or double occu- 
pancy are the same as each other between the Hubbard 
model and the mapped AM. When the virtual processes 
allowing empty and double occupancies, from which the 
RVB mechanism arises, are considered, the probability is 
0[t'^/U^'\ in the Hubbard model. On the other hand, we 
assume on the basis of Egs. p.6p . (jAlOp . and (|Alip that 



A(£) 



\t*\/ina,), \e\<W*/2 
0, |£| > W*/2 ' 



(A13) 



for the mapped AM. The probability is O[A{0)/U^] x 
0{W*) in the AM. The two probabilities should be equal 
to each other: 



0[<VC/2] - O[A{0)W*/U^ 



(A14) 



When U/\t\ ^ 1 and 1 < Wg < 2 are assumed, it follows 
from 

Eqs.lSl, (1^^ - and that 

4>^^0{U/\J\), (A15) 
it follows from Eq. dAZl and (|A15|) that 

\t*\^0{U), (A16) 
and it follows from Eqs. and (|A10I) that 

p('=)(0) = l/[7rA(0)] = 0{1/U). (A17) 

Then, 

(A18) 



lim p(0) ^ 0, 

(7/|t|->+oo 



and 



hm A{0)/\t\ 

f7/|*|->-l-oo 



+CX3. 



(A19) 



An MI transition occurs at [/c/|^| — +oo as a function 
of U/\t\. The MI transition is a discontinuous one, as 
discussed in Sec.lIVI 

Since W* = 2\t*\/^^ = 0(| J|), Tk is estimated to be 



0{\J\). 



(A20) 



If the RVB term is once self-consistently considered, even 
if (5 = is assumed, Eqs. (jAlQI) and (jA20p guarantee that 
Tk is nonzero even for U/\t\ ^ 1, though Tk K as 
U/\t\ +00. 

According to the analysis above, it is plausible that, in 
general, the GS can never be an insulator provided that 
the virtual processes allowing empty and double occu- 
pancies are possible, i.e., provided that U is finite. 

It is straightforward to show that anomalous terms 
proportional to 

e[\n{e + iO) + \n{-e-iO)] = 2£ln|e| - i7r|e|, (A21) 



appear in the multisite self-energy higher order in Js{q)- 
Since the anomalous terms are continuous at e = 0, they 
can cause no gap to open in p^'^\e) or p(e). Within a pre- 
liminary study, any anomalous term can be found such 
that it causes a gap or a zero-gap to open. Since Eq. 
()A2ip is anomalous at e = 0, the half-filled GS in one 
dimension is never a normal FL. If no gap really opens 
as discussed in Sec. lIVi it is an RVB type of TL liquid. 



Appendix B: Irrelevance of the Kondo effect in a 
continuous model 



We consider a continuous model in one dimension: 
1 f^^ 

+ dxi;t{x)ijl,{x)i;-a{x)i;a{x), (Bl) 

^ „ Jo 

where tpl{x) and ipa{x) are fermionic field operators. 
Since we are interested in the correspondence of this 
model to the Hubbard model, Hvf = 2\t\a, xi — La, 
kp — 7r/a, and gi = Ua are assumed. When the periodic 
boundary condition is assumed and 



V', 



(B2) 



is used, it follows that 



('i+qa"'l-q-aap-aaka- (B3) 

kpqa 



91 

2La 



When the continuous model (jBl|) or (jB3| is extended 
into a two-band continuous model in which the disper- 
sion relation is given by Eq. (|4.6p . it is the Tomonaga- 
Luttinger (TL) model. When U is nonzero in the Hub- 
bard model, electrons and holes with energy je — /z| 
are virtually excited in the GS; (nka) = 0(1) for k such 
that < E{k) -fi<U, and (1 - Uka) = 0(1) for k such 
that —U < E{k) — fi < 0, where nka- is the number opera- 
tor for the Hubbard model in the wave-number represen- 
tation. On the other hand, Eq. (|4.6p can only accurately 
describe electrons in the vicinity of the Fermi level, i.e., 
electrons with |fc — kp\a ^ 1 or |fc -I- k^la ^ 1. Thus, 
the correspondence is only accurate in the weak-coupling 
case of U/\t\ = 2g/{hvp) ^ 1. It is questionable how ac- 
curately and relevantly electron correlation in the strong- 
coupling Hubbard model with U /\t\ ^ 1 can be treated 
by the strong-coupling TL model with g/{hvp) 3> 1/2. 
When U/\t\ > 2 in the Hubbard model, for example, 
either of the splitting into the upper and lower Hub- 
bard bands, the superexchange interaction, and the RVB 
mechanism is a relevant effect; however, neither of them 
can be treated by the TL model with g/{hvp) ^ 1. 
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In general, when the dispersion relation is linear in k 
in a model such as one defined by Eq. (|B3I) , the model is 
essentially a continuous model. The TL model is essen- 
tially a continuous model, even if umklapp processes are 
considered as one of lattice effects. 

We denote the Green function and the self-energy in 
the continuous model (jBip or (jB3l) by Ga{isi,x — x') and 
So- (iff/, a; — x'), respectively, in the real-space represen- 
tation. When we follow the definition of the single-site 
diagram for the Hubbard model, the single-site or local 
diagram should be defined as a diagram such that it in- 
cludes only (7-lines and local lines of the Green function, 
which is simply Ga{i£i,x — x' = Q). If at least a mul- 
tisite or non-local line of Ga {isi , a; — x' ^ 0) is included 
in a diagram, the diagram is a non-local one. Thus, the 
self-energy can also be decomposed into the local Scr(ie/) 
and the non-local AYi^{iei, x — x') in a formal way: 

'i^a{i£UX~x) = Y.cr{i£l)5x-x' + AS^r (lE/ , X - ) , (B4) 

where 6x is a Kronecker-like delta defined for real x by 

It should be noted that 5x is never the delta function. 
Since a: is a continuous variable, (5a; = in a practical or 
physical sense. It is obvious that the local T,„(iei)5x-x' 
can play no role in the continuous model. The Kondo 
effect can play no role in any continuous model. 

The TL model is essentially a two-band continuous 
model. Thus, the proof of this paper for the Hubbard 
model that no complete gap can open because of the 
Kondo effect in the presence of (5 = 0+ never contradicts 
the opening of the charge gap in the TL model or the 
boson model, where the Kondo effect is irrelevant. 

Appendix C: Instability of the Mott insulator in 
infinite dimensions 

When the Hubbard model in D dimensions is studied, 
e.g., t/^/D is substituted for t in Eq. (|2TT|) . Then 

2t ° 

£;z5(k) =-^^cos(fc,a), (CI) 

is substituted for Eg. (|2.3p . The effective bandwidth of 
EnOi.) is 0(|t|) for any D. The single-site S(j(ie;) is of 
the leading order in 1/D, the multisite AI](j(i£;,k) is of 
higher order in 1/D, and the conventional Weiss mean 
field (MF), which is a multisite effect, can be of the lead- 
ing order in 1/D.— If the single-site term is rigorously 
considered and no multisite term is considered in a the- 
ory, the theory is within the S^A.^^^^^ The KLT is in a 
sense 1/D expansion theory based on the S'^ A to include 
multisite terms. 

First we consider the half-filled GS in the S^A. Since 
no conventional Weiss MF is considered, no symmetry 



can be broken in the S^A. Then, the whole analysis in 
Sec. lIIIl is also valid within the S^A. If 5 = 0+ is assumed, 
Tk > OK and the GS is a normal FL, though Tk may be 
infinitesimally low as discussed below. According to Eqs. 
(ITTTal) . (I2J51) . and (l2lB . 

k 

Thus, p{0) does not depend on U.^^ Because of Eq. p.6p . 
A(0) = l/[7rp(0)] > does not depend on U either. The 
self-energy can be expanded as in Eq. p.2p . According 
to the FL theory,— the specific-heat coefficient and the 
static homogeneous spin susceptibility are given by 

7= (2/3)7r2fc2^^p(0), (C3) 

and 

[Xs(0,q)]|q|^0-2^,p(0), (C4) 

respectively. A discontinuous MI transition occurs at 
t^c/l^l = +00 between a normal FL and the prototypic 
Mott insulator in case of 5 = O"*" . 

According to theories based on the DMFT^^— on the 
other hand, a discontinuous MI transition occurs at finite 
U in case of S = 0. When U increases, a complete gap 
opens at Uc2 = 0{\t\); both of 7 and [xs(0, q)] di- 
verge as J7 — > Uc2 — 0. When U decreases, the complete 
gap closes at Uci — 0{\t\), which is smaller than Uc2- 

In the DMFT, all the single-site terms are rigorously 
considered as the dynamical MF, and any conventional 
Weiss MF, which is a static MF, is not considered. Thus, 
the DMFT and the KLT are exactly equivalent to each 
other within the S^A. According to the study based on 
the KLT in this paper, the discontinuous MI transi- 
tion characterized by two t/c's, Ud and Uc2, must be 
one between a normal FL, which is characterized by 
Tk > K, and the Mott insulator, which is character- 
ized by Tk = K. The fact that no order parameter can 
appear in the DMFT and both of 7 and [Xs(0, q)]|q|_j.o 
diverge as C/ Uc2 — is a proof or at least a strong piece 
of evidence that Tk = K and residual entropy remains 
in the insulting GS or the Mott insulator in the DMFT. 

According to the study in Sec. lIIII of this paper, Tk > 
K for (5 > 0. Then, the following two cases are possible 
for the dependence of Tk on S: 

lim Tk > K, (C5a) 

<5-i.O 

and 

limTK = OK. (C5b) 

Based on the result of the DMFT for (5 = 0, it is antic- 
ipated that Eq. (IC5aP is satisfied for U < Uc2 while Eq. 
(jC5bp is satisfied for U > Uc2- It is desirable to confirm 
this dependence of Tk on 5 and J7 by a theory based on 
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the DMFT. li S = 0+ is assumed in the theory, the GS 
of the theory is not the Mott insulator but a normal FL. 

Since a normal FL with Tk = +0K behaves like the 
Mott insulator in a high-T phase at T > Tk = +0K, 
there is no difference between the normal FL with Tk = 
+0 K and the Mott insulator, excepting at T — OK. 
There is no practical difference between physical proper- 
ties of the phase with Tk = +0K and those of the phase 
with Tk = ok. There is no significant inconsistency be- 
tween the theoretical result based on the DMFT—"— and 
the analysis based on the KLT in this paper. 

Next we consider the half-filled GS beyond the S'^A; 
but only the RVB term is considered beyond it within 
the KLT, as in AppendixEl When U/\t\ > 1, the 
superexchange-interaction constant is given by 

Jd = -^tyiDU), (C6) 

between nearest neighbors. The stabilization energy by 
the RVB mechanism is 0(| J/jj) per unit cell, so that 

kBTK^O{\JD\)^0[ty(DU)]. (C7) 

The RVB mechanism is of higher order in 1 /D. Following 
ApDendixf^ it is straightforward to show that p(0)|t| — >■ 
as U/\t\ — > +00 for any finite D; p{0) is not given by 
Eq.JCSl). If P(.0)\t\ < 1 for finite D and U/\t\ > 1, the 
RVB mechanism is crucial, so that the GS is simply an 
RVB type of normal FL. On the other hand, if the limit of 
D +00 is followed by that of ;7/|i| ^ +oo, the RVB 
term vanishes, so that p(0) is given by Ea. (IC2p . The 
constancy of p(0) as a function of U is not an property 
o{ D —i' +00 dimensions but one in the S'^A. 

Thirdly we consider the true half- filled GS. When 
U/\t\ ^ 1, the Neel temperature Tn is given by 

kBTj, ^ D\Jd\/2 = 2tyU, (C8) 

in the MF approximation, which is rigorous in the limit 
of D — 7- +00. The Weiss MF of magnetism is a leading- 
order effect in Since 

Tn »Tk -0[iVPt^)], (C9) 



the GS is a local-moment type of AF magnet. 

In a high-T phase at T ^ Tk, the susceptibility of the 
mapped AM obeys the Curie- Weiss (CW) law: 

X.(0)~l/[fcB(T + TK)]. (CIO) 

According to the KLT, the effective Fermi energy Ep 
discussed in Sec.|T]is as large as Ep ~ /cbTk. In a high-T 
phase at T > Tn and T > E^/kB or T > Tk, 

X.(0,q) = X.(0)/[1 - (l/4)J.(q)x40)] 

~ (l/fcB)/[T + TK- J.(q)/(4fcB)], (Cll) 

where 

D 

J,(q) 2Jd J2 cos(g^a). (C12) 

1^=1 

The susceptibility obeys the CW law for any q because 
of the T dependence of Xs(0); the Weiss constant de- 
pends on q. Since electrons behave like local moments, 
the high-T phase can be regarded as the Mott insulator. 
When U/\t\ > 1, the GS is an AF type of Mott insulator. 
When U/\t\ ^ 1, the GS is also an AF insulator. Since 

Tn « Tk = 0(K|/fcB), (C13) 

it is an itinerant-electron type of AF magnet. The spin 
susceptibility obeys the CW law at Tn < T <C Tk for 
particular q such as q ~ (7r/a)(±l, ±1, • • • ,±1) because 
of the nesting of the FS. When U/\t\ < 1, the GS is not 
any type of Mott insulator. 

Lastly we consider the half-filled or non-half filled GS 
for U/\t\ > 1 and finite D within a treatment where the 
single-site T,a{i£i) is rigorously considered, the RVB term 
is considered, all the other multisite terms are not con- 
sidered, and no Weiss MF is considered. Within this 
treatment, the GS is an RVB type of normal FL. This 
type of normal FL is so useful that it can be used as a 
unperturbed state to study high- Tc superconductivity in 
the vicinity of an AF type of Mott insulator 1^ 
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